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A Fourier analysis of parametric level dynamics for random matrices periodically depending on 
a phase is developed. We demonstrate both theoretically and numerically that under very general 
conditions the correlation C{ip) of level velocities is singular at = for any symmetry class; the 
singularity is revealed by algebraic tails in Fourier transforms, and is milder, the stronger the level 
repulsion in the chosen ensemble. The singularity is strictly connected with the divergence of the 
2nd moments of level derivatives of appropriate order, and its type is specified to leading terms 
for Gaussian, stationary ensembles of GOE, GUE, GSE types, and for the Gaussian ensemble of 
Periodic Banded Random Matrices, in which a breaking of symmetry occurs. In the latter case, we 
examine the behaviour of correlations in the diffusive regime and in the localized one as well, finding 
a singularity like that of pure GUE cases. In all the considered ensembles we study the statistics of 
the Fourier coefficients of eigenvalues, which are Gaussian distributed for low harmonics, but not for 
high ones, and the distribution of kinetic energies. 
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I. INTRODUCTION 



In this paper we shall be concerned with spectral properties of random matrices H{ip), periodically depending on 
a phase (p ejO, 27r]. Random matrices j^] are often used to describe statistical properties of disordered mesoscopic 
systems [^HJj]. In particular, random matrices parametrically dependent on a phase are used as models for conduc- 
tion in small metallic rings threaded by an Aharonov-Bohm flux wh/e, and great attention is being devoted to the 
dependence of their eigenvalues ei{i.p) on the phase ip |||,^,0,||J|,|l3l ■ Objects of particular interest are the statistics 
of level slopes e^(<^), (also called level velocities; a prime denotes a derivative with respect to f), of level curvatures 
e'-{ip), and velocity- velocity autocorrelations 

Civ) ■■-^J + v)) V = ^ , (1-1) 



where ( ) denotes ensemble averaging, the bar denotes phase averaging, and A is the average level spacing. 

Autocorrelation functions that are closely related to (LI) (at least insofar as a single-particle description is valid) 
are also accessible to laboratory experiments on quantum dots [ pd] | . 

To the extent that random matrices can be taken as models for the Hamiltonians of quantum systems that become 
chaotic in the classical limit, the identification of "universal" properties of parametric level dynamics is also relevant 
for quantum chaology [ p^ . In that case the phase i/j represents any external parameter which determines the spectrum. 
The motion of levels as functions of may be thought of as a dynamics of fictitious particles (with positions e^), with 
</? playing the role of "time" . 

In a series of papers, Altshuler and coworkers have shown that for disordered metallic systems the parametric 
statistics are in fact universal, after rescaling the levels to unit mean spacing and also rescaling the parameter ip as 
X = -KLp^/c, with c = C(0) denoting the mean squared velocity of energy levels. Therefore, the only two system 
-specific parameters are A and c; the latter has the meaning of twice the average kinetic energy of the level moving 
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in the time Lp. An example of such a parametric universality has been discussed earlier in studies of level curvatures 

@Oji|. 

In particular, the theory |j] leads to a remarkable scaling property of the velocity- velocity autocorrelation (1.1): 



C(V3) « C^{-KLp^c) = C$(X). (1.2) 



In Eq.(1.2) $(.) is an universal function, only dependent on the general symmetry properties of the chosen matrix 
ensembles. Determining the shape of this universal function is manifestly an important task. At large X , the behaviour 
<I>(Ar) ~ — was theoretically predicted and numerically tested instead, the small X behaviour appear 

to critically depend on the chosen ensemble. Explicit expressions have been obtained for a closely related [but 
distinct from C(</j)] autocorrelation function at fixed energy. The resulting expressions are valid for homogeneous 
ensembles only, in which the symmetry class of the random matrix is independent of (/s, and are quite complicated (for 
the orthogonal ensemble a triple integral with regularization is required). A global approximation for ^{X) has been 
proposed |[l^ ] . For the case of a classically chaotic system subject to a Aharonov-Bohm flux Berry and Keating |l^] 
have recently obtained a semiclassical approximation for C {(p) having the form of an everywhere analytic function of 



In this paper we shall investigate correlations (1.1) in various ensembles of phase-dependent random matrices, 
using Fourier analysis, which is particularly well suited to this purpose, because of periodicity, and also because 
the singularities of C((/5) are easily detected by this technique. We shall consider the case of homogeneous Gaussian 
ensembles, with orthogonal (GOE), unitary (GUE) and symplectic (GSE) symmetry, and also the physically important 
situation, in which the time - reversal symmetry is broken and a transition from GOE to GUE occurs while switching 
on the "Aharonov-Bohm flux" if. The latter case will be analyzed on the ensemble of Periodic Band Random Matrices 
(PBRM) which exhibits the three characteristic regimes of conduction in disordered solids, namely, the localized, 
the (proper) diffusive, and the ballistic one. 

A general theoretical approach to the Fourier analysis of spectral correlations is developed in Section II, where it 
is shown that, under very general conditions, C((^) must have a singularity at = 0, signalled by an algebraic tail 
in its Fourier transform. This result is valid even in cases (such as homogeneous GUE, GSE) which were hitherto 
believed to be singularity-free. The singularity is however milder, the stronger the level repulsion in the considered 
ensemble. We determine its exact type for various symmetry classes. For the homogeneous cases our predictions are 
confirmed by numerical results, presented in Section III, and by the exact solution for matrices of rank 2, presented 
in the Appendix. The model of Periodic Band Random Matrices corresponding to the Aharonov - Bohm case is 
analyzed in Section IV. We present theoretical arguments and numerical results supporting the thesis that for this 
case the behaviour of the correlation function at </? = is similar to the one observed in the homogeneous GUE case, 
and we compare our data to other theoretical predictions. Furthermore we present distributions of avoided crossings 
for PBRM, which are deeply related to the properties of the velocity correlation function, and briefly discuss the 
modifications of the correlation function which occur on moving towards the localized regime. 

Further applications of Fourier analysis presented in this paper include the statistics of various quantities. Among 
these, the Fourier amplitudes of level velocities (currents), which are found to be Gaussian for low harmonics but not 
for high ones, and the mean squared current, which appears to follow a distribution with an appropriate number 
of degrees of freedom. Finally, we discuss the moments of level derivatives of order k with respect to phase and 
find them either to diverge (if k is not small enough) or to be directly related, via the scaling law, to the average 
kinetic energy of levels. 



II. FOURIER ANALYSIS. 



In full generality, by a Periodic Random Matrix (PRM) we mean a random Matrix H{ip) periodically depending 
on a phase ip; more precisely, a matrix- valued, periodic, stochastic process, taking values in a class of matrices having 
a definite symmetry (symmetric, complex- hermitean, symplectic), which is not necessarily the same for all values of 
Lp. An ensemble of PRM is not an ensemble of matrices, but an ensemble of trajectories H{lp) in a space of matrices; 
every such trajectory will be called a realization of the ensemble. In this paper we shall consider stationary ensembles, 
such that their statistics is invariant under translations mod(2n) of the phase ip, and non-stationary ones as well. 
Although the numerical results described in this paper were obtained for some very specific, well-defined ensembles, 
the argument to be presently described is applicable under rather general assumptions, hence to a much broader class 
of ensembles. 

Our first assumption is that, with probability one, the trajectory H{ip) is an entire function of ip. On the strength of 
this assumption, the eigenvalues e{ip) of almost any realization H(ip) are analytic functions of ip, and their behaviour 
in a neighborhood of any point ip in [0, 2tt] can be described by analytic perturbation theory. The corresponding 
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convergence radii, however, are reahzation-dependent. If they can become arbitrarily small in a statistically significant 
fraction of the ensemble, perturbative methods will not be applicable to the calculation of certain statistical averages. 
This is precisely the situation we are going to discuss below. 

Given a realization H(ip), any eigenvalue e{ip) can be expanded in a Fourier series. The Fourier expansions we shall 
use are : 

+ 00 +C30 

e(^) = ^(^) = E c^e"'"' (2-1) 



where a„ = al„ and, on account of Eq.(Ll), c„ = n A \\o.n\ y . The statistical average (.) is meant as an average 

over all the eigenvalues lying in a selected energy range (in the following we drop a subscript i denoting a given level 
in Ci), followed by an average over all r ealizations in some ensemble. Because C{ip) is a correlation function, it must 



be even, C{ip) = C{2tt — tp). From (|2.lD it follows that 



+00 



= C(0) = ^c„. 



If, in the chosen ensemble, a scaling law of the form ( |l.2| ) holds, then, at c >> 1 : 

Tl 

Cn « V^P(^), (2.2) 
with P{t) a function related to the Fourier transform of <i>(X) [Eq.( |l.2| )]: 

P{t) = (27r2)-i / $(X)e-'*^/'^dX (2.3) 



It follows that P{t) defines a probability distribution on the real line, with (not necessarily finite) moments A^°''> = 
In particular, the width of the Fourier spectrum is measured by: 

+cx> 

E \n\cn 

^ ^Vcj\t\ P{t)dt = (2.4) 

— 00 

Even moments are related in an obvious way to derivatives of $ at the origin: 

$(2'c)(o) = (_i)fc7r-2fc^(2fe) (2.5) 

In ensembles for which the scaling theory is valid, the functions (f>, P and the constants ^^"^ depend on the 
universality class only, which will be occasionally specified in the following by a suffix (3 (as, e.g., in $^), with /3 = 1,2,4 
for the GOE, GUE, GSE universality classes, respectively; or by the suffix AB in the symmetry-breaking, GOE-GUE 



case. Eqs.(p^),(2.4) provide a way to check the scaling hypothesis in the Fourier domain. For the specific ensembles 
considered in the following Sections they were indeed confirmed by our numerical computations, to be discussed later. 

Following our previous assumption, H{ip) can be analytically continued to the whole complex t/?— plane. Any eigen- 
value e{Lp) is an analytic function of the complex variable tp, except for branch-points (BP) Zj = ipQj+ijj corresponding 
to complex level crossings The index j labels the BPs, j — 1, ...TV, with M the total number of BPs exhibited by 
the given level as a function of (p in the complex domain. Since e{ip) is real at real (^pthese BPs are distributed sym- 
metrically with respect to the real axis, but none of them falls on the real axis itself El. Therefore, given a realization 
H{(p), the function e{ip) is analytic inside a strip of halfwidth P around the real axis, P being the smallest of the \ 
In order to estimate the complex Fourier coefficients 



^Although crossing of different eigenvalues at real values of ip is possible, such real crossings would not spoil the analiticity of 
e{ip). 
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with n > we use the multiple path illustrated in Fig.^ In the following we assume n > (if n < the path 
must be reflected with respect to the real axis; the corresponding analysis differs from the following one by obvious 
modifications). In this way 



where Ij is the complex contribution of the j — th branch-line: 



(2.6) 



CXD 

-inzj 



Ij = J dx e-"-^e,(x), (2.7) 



and Sej{ x) is the jump of e{ip) across the j— th branch cut, at the point ip ^ Zj + ix. The behaviour of the Laplace 



integral (2.7) at large n is determined by the behaviour of 5ej{x) close to the BP, as we shall presently discuss. We 
concentrate on BPs with 7^ << 1, whose contribution will dominate at large n. We shall consider simple BPs only; 
in other words, we shall attribute zero statistical weight to crossings involving more than two levels. If a realization 
H{ip) (if complex) is thought of as a two-(real) parameters family of matrices, then generically such a family will not 
contain matrices exhibiting higher-order degeneracies, because multiple level crossings cannot generically be achieved 
by varying just two real parameters. Therefore, the assumption we need to justify our neglecting higher-order BPs is 
just that such non-generic realizations H{ip) have zero probability; this is a smoothness assumption on the statistical 
distribution of the ensemble, which appears fairly general. 

For ip real and close to ipoj the considered level is undergoing a narrow avoided crossing. It is in fact close to 
degenerate with a different level, the two eigenvalues being approximately: 

e± (if) « Ej ± h, y'7| + (^-^o,)2. (2.8) 



Carrying this to complex ip we get, for cp close to Zj, 6ej{(p) ~ 2hj^jjX. Substituting into (2.7), and using known 
estimates for Laplace integrals [Q, we obtain 

^k/^n'^^^e-'"'^. (2.9) 

TT 

We now introduce the single-level velocity correlation function (CF): 

+00 

R{^) :- A-2e'((^')e'(¥'' + ^) = ^"S"'"' (2.10) 



(the bar denoting average over tp'), which depends on the chosen eigenvalue, and on the realization as well. This 



fimction is connected to the original CF by C{ip) — {R{ip)), and c„ = (r„) . From (2.9) we obtain: 



n 

Thus, prior to ensemble avera ging R(ip) has an exponentially decaying Fourier expansion, and is therefore analytic in 



a strip of halfwidth 2r. From ( 2.11 ) it actually follows that R{lp) has logarithmic singularities ai ip — ±2ir-f 2fc7r, with 
leading terms F log((^±2ir-|-2/c7r). Therefore, for real (p close to 0, R{(p) behaves like F log(</j^-f4F^). On averaging over 
the matrix ensemble, arbitrarily small F come into play, that is, the logarithmic singularities of R{p) come arbitrarily 
close to the real axis. As a result, C{ip) will not be analytic at = any more; the type of singularity will depend 
on the statistical weight which is given to small F in the chosen ensemble. Actually, a standard probability-theoretic 
argument, quite independent of the above elaborations, shows that in order for C{ip) to be analytic, too small F must 
not be allowed at all! 
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In fact, if we assume C{(p) to be analytic, then, being also periodic, it must have an analyticity strip of (strictly) 
nonzero width 2L > 0. Then the coefficients c„ must satisfy 

for all L' < L (with appropriate constants Kl/). Since c„ — 7i^A^^(|a„p), the probability pn{L") that |a„| be larger 
than I"' will satisfy the Chebyshev inequality pO[ |: 

p„(L") < c„n-2A2e2L"|„l < if^,A2e-2(i'-i")l"l. 

Therefore J2nPn{L") < +00 as soon as L" < L. Then the Borcl-Cantelli Lemma of Probability theory ensures 
that, with probability one, all but finitely many of the coefficients a„ satisfy |a„| < I"' for any L" < L. In turn, 
this implies that e{(p) is analytic in a strip of width L with probability one; that is, the probability of a branch-point 
lying closer than L to the real axis is zero. Therefore, in order to get an everywhere analytic C{ip), the ensemble must 
be such that the distribution of the imaginary part of BPs has a finite gap. Since the latter distribution is clearly 
related to the distribution of the widths of avoided crossings, this appears to be rather an exceptional situation. 

In order to get more precise information about the nature of the singularities ofC{(p), a knowledge of the statistical 
distribution of the BPs is needed, which cannot be obtained under the general assumptions made up to now. We shall 
then assume that this distribution is described by a density p{'~f,(po), and also that different BPs give uncorrelated 
contributions. Furthermore, in the rest of this Section we shall restrict ourselves to stationary ensembles, such that 
all statistical properties are independent of (p; explicit examples will be given in the next Section. In such cases 
we can assume statistical independence of (po and 7, and a uniform distribution for the former. In other words, 
p(7, ip„) = /(7)/(27r). Then, from Eqs.(|2j) and (|2j) we obtain 

- / 7e-'"V(7)- (2.12) 



In order to assess how does this expression scale with c, we first notice from (2.8) that the quantity h sets the scale 
for the level velocity close to the avoided crossing. Therefore we assume {^hyA^'j to scale as the variance of the 
rescaled level velocity, that is, as c = C(0) = ^c„. At the same time, /17 is the width of the avoided crossing, and 
this identifies the scale of 7 with that of A/h ~ c~^/^. On such grounds, setting 7 — xc^^^^, we expect 



/(7)=cl/2^^(^cV2), 



with F(x) a scale-independent distribution. Substituting this into (2.12) we obtain 



dxF[x)xe 



-2x7i/ \/c 



(2.13) 



(in the rest of this Section, a will denote undetermined numerical factors). In order to estimate A/", we identify its 
order to magnitude with that of the number of avoided crossings experienced byefw) as varies between and 27r, 
which has been found to be proportional to ^Jc using two different approaches |22|,|24|. We are thus finally led to 



00 

a— / dxFix)xe~ 
n J 



■2xn I -^/c 



(2.14) 



This formula has a number of implications. In the first place, it is exactly of the form (2.2); the above argument has 
led us to an independent confirmation of the scaling law at n >> ^/c. Second, it yields an asymptotic estimate for 
the universal function P{t) at large t : 



00 

P(t) ^ J J dxF{x)xe'' 



(2.15) 



It does not appear legitimate, ho wev er, to use formula (2.14) outside the asymptotic regime n >> ^/c, because of the 
underlying asymptotic estimate (p^). 
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To pr oceed further, we need to specify F{x), or at least its behaviour at a; << 1, which determines the asymptotics 
of ( ^11 at large n. The variable x is, apart from rescaling, the same as 7, which is in turn closely related to the width 
of an avoided crossing. For the case of Gaussian ensembles, the statistics of the latter quantity has been studied in 
p2| , ^j24| , where the probability of small avoided crossings of size e was found to scale as e^~^, with /3 the repulsion 
parameter characteristic of the given ensemble. This result is actually very general. According to (2.8), the probability 
of a small level spacing less than e, which scales as e^^^, is also proportional to the probability of getting a BP inside 
a small circle of radius e in the complex 93— plane. The latter probability scales as e Jq fdj, whence /(e) ~ e^~^. 

Our last assumption will be that level repulsion in our ( stationary) ensemble is ruled by a repulsion exponent /3. 
Then, substituting F{x) 



( at small x) into Eq.(2.15|) we finally get 



P{t) - at 



-2-13 



(2.16) 



which is the main result of this Section. It implies that in no stationary ensemble can Ci3{ip) be a smooth function 
at if — 0, because its (/3 + l)-th derivative must have a singularity there. It sho uld b e noted that the asymptotic 
behaviour c„ ~ n~'^~^, which leads to this conclusion, directly follows from eqn.( ^.l^ ) and fr om t he fact that the 



probability of small 7 scales as 7^"^, independently of the sca ling c onsiderations following eqn.(2.12). 

The nature of the singularity associated with the behaviour ( ^.16 ) now follows, either from known results in Fourier 
analysis, or, more directly, from statistically averaging the previously described small-(y9 behaviour of R{(p). The 
resulting small-y? expansions for the cases P — 1, 2,4, written in terms of the scaled variable X = nipy^, are (up to 
leading terms only): 



Ci((p)/c - l + 6iX2ln|X|, 

C2{'p)/c^l~b2X^ +b3\X^ 

Ci{ip)lc - 1 - biX^ + 65^4 + 66 



(2.17) 



Some of the coefhcients 6„ will be specified later, for specific choices of the PRM ensemble. 



A. Scaling and Level Derivatives. 



The statistics of BPs not only determines the statistics of Fourier coefficients, but also those of the curvature 
K{ip) = e"{ip) and of higher-order derivatives as well. Several results are known about the distribution of curvature in 
the case of stationary ensembles [^3|,^ 25 2^. Here we shall point out some immediate conclusions that can be drawn 
from Fourier analysis, concerning the statistics of an arbitrary lev el derivative Ka{'~p) = e''°'\ip) (being analytic, levels 
have derivatives of all orders), under the scaling assumption (1.2). 

From the very definition we get the mean square value of such a derivative in the form: 



+00 



(2.18) 



where, as usual, the bar denotes an average over the phase (p. On ensemble-averaging this equation we get 

^ A2^n2("-i)c„ (2.19) 



From the scaling law we now obtain: 



+00 



(2.20) 



where At^^-^) is the {2a~ 2)-th moment of the distribution P(t). From this equation, using the relation (2.5) between 
the moments of P{t) and the derivatives of C{lp) at ip = Q, we obtain that the (2q; — 2)-th derivative of C{(p) is singular 
at if = if, and only if, the 2nd moment of the a-th derivative has nonintegrable singularities as a function of if. 
On account of ( |2.17D , this occurs when a is larger than 1,2,3, respectively, for GOE, GUE, GSE-like level repulsion. 
If, in addition, the PRM ensemble is stationary, then, fo r suc h values of a, the moments are divergent at all ip. For 
a = 2, that is, for ordinary curvature, and for (3 > 2, Eq.( 2.20|) establishes the proportionality of rms curvature to the 
average kinetic energy: 
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where (^.5|),(2.17) have been used to find A^'^\ This relation is similar to one found by Akkermans and Montambaux 
p?! for the case of a symmetry-breaking AB flux, the difference being the additional (^-average performed here. If, in 
addition, the ensemble is stationary, the same relation holds for the rms curvature at any point (p. 

III. STATIONARY ENSEMBLES 

A. Spectral correlations 

The stationary ensembles we consider are defined by: 

H{(p) ~ Hi cos tp + H2 sin (p, (3-1) 
where Hi,H2 are fixed, (i.e., (/?— independent) random matrices, drawn independently from one Canonical Gaussian 



ensemble (GOE, GUE, GSE). Depending on the chosen ensemble, the PRM ensemble defined by (3.1) has trajectories 
in different matrix spaces (the space of real symmetric matrices of rank N for GOE; of complex Hermitean matrices 
of rank N for GUE; of symplectic matrices of rank 2N for GSE). For every tp, the elements of the matrix H{p) are 
centered Gaussian random variables, of variance 

'(i?i,2).,f) = ^(l + %). (3.2) 

With the chosen normalization of matrix elements, the average eigenvalue density of H{ip) obeys (in the N ^ 00 
limit) the Wigner semicircle law in the form 

p{e)^^{4-e'f\ (3.3) 

In the central region of the spectrum the average level spacing is then A = p~^(0) — ir/N. 



Because the equation (3.1) resembles the parametric equation of an ellipse, the PRM ensembles just defined will be 
termed "Gaussian Elliptic" in the following. They have a number of distinguished features. In the first place, they 
are stationary ]l^ , |l5t ; therefore, all "pointwise" distributions, such as the distribution of velocity e'{(p), of curvature 
e"{ip), at given ip, are independent of (p. Moreover, two-points statistical correlations such as {e'{ipi)e'{ip2)) only 
depend on the difference \(p2 — ipi \ ■ 

Secondly, at any ip the derivative H'(ip) is statistically independent of H{ip), and belongs to the same canonical 
ensemble; indeed, a simple computation shows that the matrix elements of both H{ip) and of H'{p) are Gaussian 
distributed, and that they are uncorrelated. As a consequence, the level velocity e'{(p) has a Gaussian distribution. In 
fact, e'{ip) is given by an expectation value {e\H' {ip)\e) on the eigenvector |e) of H{ip) corresponding to the eigenvalue 
e{ip). On account of the independence of H{p) and H'{ip), and of the rotational invariance of the Gaussian ensembles, 
such a matrix element must have the same distribution as any diagonal element of H'{ip). Therefore it has a Gaussian 
distribution, with variance 

An independent (albeit more complicated) pro of o f the Gaussian distribution of level velocities can be obtained from 
supersymmetric calculations 0] . The variance ( 3^ ) also follows from a statistical mechanical analysis of level dynamics 



(Eq.(3.8) of [p^) u nder the assumption of the applicability of the canonical ensemble. 

Combining ( |3.4D with the mean spacing A we obtain for the center of the spectrum the average variance of level 
velocity (also called the average squared current) 



2 



2^. (3.5) 



We have numerically computed Fourier transforms of eigenvalues e((p), for matrices of the above type, of rank 
TV = 40 ^ 400. The eigenvalues were unfolded in order to normalize them to unit level spacing. In computing (fast) 
Fourier transforms we have used grids of 128 ~ 1024 points in (0, 2tt) - the optimal size of a grid grows with the 
matrix rank N. To avoid confusions, in the following the results of averages over (p taken for a fixed level and a fixed 
realization will be called " mean values" , while by " average values" we shall mean results of statistical averaging (over 
realizations, and/or over different levels). Thus "mean" values will be fluctuating quantities, in contrast to "average" 
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valu es. F or each level we have computed Fourier coefRcients a„, mean correlation Fourier coefRcients r„ [defined by 
Eq.( ^.ll" )], and the mean squared velocity r = ^rn. Finally, statistical averages were taken over all the eigenvalues 



lying in the central half of the spectrum, and the results were further averaged over samples of 20 300 matrices. 



As a first step, we have checked the scaling law (O). The numerically computed moment defined by (p.4|), 

grows as VN, hence, predicted by the scaling law. Moreover, the relative width of the Fourier spectrum, 

measured by 5n := {{n^) — {\n\)^y^^ / {\n\) is constant, as required. Direct evidence for the scaling law in the GUE 
case is presented in Fig.0, where the function P2{t) is plotted for various matrix sizes N against the rescaled variable 
t = n/{\n\). 

All the assumptions which underlie the argument developed in the previous Section app ear t o be satisfied for Elliptic 
ensembles, so that the small- ip behaviour ( 2.17| ) or, equivalently, the large-i behavio ur (2.16| ) are expected. 



Some of the coefficients 6, which belong to the nonsingular part of the expansion ( 2.17] ), are explicitly known 



62 = 2, &4 = 3 (3.6) 

as can be obtained pq | by expanding C{(p) around (p = to the lowest nonvanishing order, and using for the variance of 
curvatures, {e" {(p)'^^the values derived from exact curvature distributions ||r^ , |25| , p6|| . The coefficient 62 has also been 
found via supersymmetric calculations 0|. In the case of rank-2 matrices we have exactly and explicitly calculated 
the functions Ci, C2, which are reported in the Appendix and have the predicted small-iy9 behaviour. 

In numerically analyzing the large-n behavior of the coefficients c„, two competing numerical effects, which obscure 
the real asymptotics, have to be taken into account. The first is "aliasing", an artifact originated by the discretization 
used in Fourier transforms, which tends to bend the numerical tails upward. A second effect, which tends to bend 
them downward, is due to finite statistics: however large the sample of matrices, the avoided crossings sampled in 
it have a finite nonzero minimum width, that will be detected by the Fourier transform as soon as n becomes large 
enough. To get rid of such effects, we have just restricted attention to that part of the c„ sequence that proved stable 
against increasing the Fourier transform basis and the size of the sample. Results are shown in Fig.H for the three 
cases (3 — 1,2, 4, and agree with the prediction c„ ^ n^"^^^ . 

Finally, we have used the numerical Fourier transforms to reconstruct the scaling functions <I>^(X). Results are 
shown in Fig.^ for the three Elliptic ensembles, and in more detail in Fig.| for GUE . In this case we have compared 



the empirical $2 with the theoretical predictions, eqs.( 2.17 ) at small X, and X ^ aX large X, finding a good agreement 
in both cases. 



B. Statistics of Fourier coefficients 



An arbitrary rotation ip \ — > p + a results in multiplication of the Fourier coefficients a„ by the phase factor e™". On 
the other hand, such a rotation must leave the statistics of the stationary ensemble unchanged. Hence (a„) = e*""(a„), 
and, since a is arbitrary, (a„) = 0. A quite similar argument shows that {a*firn) = if n ^ m, which means that 
the complex variables a„ are pairwise uncorrelated. Nevertheless, they cannot be statistically independent. In fact, 
for any fixed realization they decay exponentially fast at large n, but their variances decrease algebraically instead. 
This sharp contrast between the individual and the average decay entails a strong statistical dependence between 
coefficients Qk a t large n. 

Formula ( ^.6[ ), shows that are sums of a large number \fc ^ \fN of contributions Ij, that come from the 
different BPs. On assuming these contributions to be independent, one would naturally expect a Gaussian distribution 
for a„, which implies an exponential distribution {x^) for |a„p. Numerical computations have confirmed this pre diction 
for all three universality classes and n = 1, 2, 3 with N = 50. However, at very large n the asymptotics ( 2.6,2T9|) holds, 
with only a few terms surviving in the sum, corresponding to very small, and unlikely, avoided crossings. In this case 
a Gaussian distribution can hardly be expected. A Gaussian distribution of all the Fourier coefficients can be in fact 
ruled out: for otherwise these coefficients, being uncorrelated, would also be independent, a case that we have just 
excluded. In addition, other variables linearly depending on the a„ such as, e.g. the curvature K{(p) = — ^ 7T,^o„e*"''' 
should have a Gaussian distributioncl, which is known to be false 1 15|, |25| , p6| | . 

Our numerical data have revealed broad, non Gaussian tails of the distributions of the real and the imaginary parts 
of a„ for sufficiently large n. In Fig.O wc illustrate the transition from Gaussian to non-Gaussian distribution of a„. 



^In fact, in the GUE, GSE cases the asymptotic behaviour of {|an|^) entails that the series for K is mean-square convergent, 
so that a Gaussian distribution of all terms in the series would enforce a Gaussian distribution of its sum. 



8 



which occurs as n increases, by plotting the logarithmic variance (i.e., the variance of ln(|a^|)) versus t = n/(|n|),for 
the case of Elliptic GUE. The logarithmic variance corresponding to a Gaussian distribution of a„ is equal to 7r^/6; 
this value is indeed found at small n < (for N = 40, w 12), but not at larger n, yielding evidence that the 
distribution of Fourier coefficients is not Gaussian any more. Moreover, the same FigJ^ shows that the log-variance 
docs not scale with n/{\n\). Therefore, although the 2nd moments of the Fourier coefficients scale in a simple way 
(1.2), uniformly (that is, for all n), the same is not true for their global statistical distribution, except perhaps in 
the range n < y^, where, due to the Gaussian character of the distribution, the scaling of 2nd moments entails the 
scaling of all moments. 



C. Curvatures, and mean kinetic energy 



Eqn.(|2!20|), and its consequences, are obviously valid in the case of Gaussian Elliptic ensembles. In addition, 
= {K^{ip)) for all (f because of the stationarity of the ensemble. From Eqs. ( ^.17 )( |2!5| )( |3.6|) one finds A^^' = 47r2 



and A^^' = 167r^/3. Thus, in particular, in the GUE Elliptic ensemble the root mean square (rms) of the curvature 
is, at any point ip, proportional to c, with a factor 27rA. 

The Gaussian distribution of the level velocities entails a Porter-Thomas distribution of squared velocities. Never- 
theless, such a distribution cannot be expected for the mean squa re level velocity r = ^ r„ (the square level velocity, 
averaged over ip; a realization- dependent quantity, see Eq.( ^.ic| )), which is not the square of a Gaussian variable. 
Braun and Montambaux [ p^ have numerically studied the distribution of r for the 3D Anderson model in the metallic 
regime, and have approximated their data by a log-normal distribution. In this Section we describe our results for 
the distribution of r in the Elliptic GUE ensemble. 

Our data show that the variance of r is proportional to N^^^, while the logarithmic variance a^^^ = (In^(r)) — 
(ln(r))^ scales as 7V~^/^. Thus the relative fluctuation of r vanishes in the limit c oo, that is, r/c is a self-averaging 
quantity. Figure displays the distribution P(ln(r)) of the logarithm of r obtained from 200 GUE matrices with 
N = 50. The small but significant asymmetry with respect to the mean advises against a Gaussian fit; as a matter 
of fact, a much better fit for the empirical distribution of r was obtained by means of a xti''') distribution, with 
the average (r) = c, and with the number of degrees of freedoms ly « '^/'^rin- latter estimate follows from the 
observation that the logarithmic variance ctJ^u of the xt distribution does not depend on the mean and is given by 
\I''(j^/2) where the Trigamma function 5''(i//2) ~ at large i/. The resulting distribution (solid line in the 
figure) fits the data much better than the log-normal distribution. A similar distribution describes the mean kinetic 
energy statistics for GOE and GSE. 

As already mentioned, afi^ depends on the matrix size as N~^/'^ . Hence v is proportional to \/N and also to the 
mean coefficient (n); in other words, the number of degrees of freedom is proportional to the width of the Fourier 
spectrum. Numerical data give for GUE v k, 3(|n|). Therefore, our results appear to indicate that the statistics of 
r = ^n^|a„p are determined only by the Fourier amplitudes a„ with n < a/c, which, as we know, are Gaussian and 
uncorrelated. This would also explain the observed dependence of the variance of r, for 

1 

-'r- E Cl^c'/'jp^m (3.7) 
|n|<yc 

This interpretation is far from complete, though. In fact the scaling law implies that the relative contribution to r 
of the region n > y/c is independent of c, so that there is no "a priori" reason why it should be negligible. 



IV. PERIODIC BAND RANDOM MATRICES 



A. General Properties 



Let us consider a circular array of N sites on a ring, labelled by the index j — 1,...N. By definition, a Periodic 
Band Random Matrix has nonzero matrix elements Hjki^p) only for sites which are no more than b sites apart on the 
circle, where 6 < A'^ is an integer specifying the width of the band. The matrix elements are further specified by: 

Hjki^p) = hjk if |j - fc| < 
Hjkif) = hjkc"^ \i j <k + h- N 
H,k{v) = h,ke~'^ ii]>k + N-h, 
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where hjk are real Gaussian variables, independent of ip, hjk = hkj, with the normalization 



(4.1) 

For an odd value of N and b = {N + l)/2 we get at ip ~ & full GOE matrix. 

PBRMs are obtained from a Bloch decomposition of infinite, periodic, banded random matrices, and tp is the corre- 
sponding Bloch index. The PBRM ensemble differs from the periodic random matrix ensembles matrices considered 
above in two important respects: 

1. H{(p) exhibits a proper diffusive regime in the range 1<<6<<7V<<6^ (which should not be identified with 
the ballistic regime b « N), and a localized regime at iV > All statistical properties of H{(p) depend on b 
and through the scaling variable x = b'^ /N, which is proportional to the ratio between the localization length 
of the eigenfunctions in the limit N = oo and the " sample size" N ; 

2. The self-adjoint matrix H{(p) satisfies H{—(p) = H*{(p). Therefore at (/? = and ip= tt the matrix H(ip) is real 
symmetric, while at any other value of ip it is complex-Hermitean. This is exactly the behaviour expected for 
an Hamiltonian depending on an Aharonov-Bohm magnetic flux. 

On account of both the properties mentioned above, PBRMs provide a model for the dynamics of electrons in small 
disordered rings threaded by a magnetic flux tp (in dimensionless units). Previous studies have shown that the 
behaviour of the ensemble-averaged, zero-flux scaled absolute curvature IC := (|e"(0)|) is quite similar to the one 
expected of the average conductance in the scaling theory of localization. The distribution of K, also depends on x, 
being identical to the generalized Gauchy distribution (first conjectured by Zakrzewski and Delande |15|| , and then 
proven for GUE and for a broad class of ensembles of random matrices |^) in the delocalized limit a; ^ oo, 
and resembling a log-normal distribution in the opposite localized limit. In general, all statistical properties change 
continuously as the localization parameter varies from the localized to the delocalized limit. For instance, Fig.^ shows 
the X— dependence of the average squared current c, the variance tr^, the log- variance cr^i^, and the moment of 
the Fourier spectrum. The latter is close to unity in the localized regime; but in the diffusive regime it scales as ~ x^^^ . 
This result has a simple explanation: the Thouless conductance (|i4r(0)|) obeys Ohm's law in the diffusive regime 
[ ||, hence it is proportional to x. On the other hand, according to a relation discussed by Akkermans and Montambaux 
[ 27t and numerically verified for PBRM |^ and 3D Anderson model [H, in that regime the Thouless conductance is 
also proportional to Cab(O). Therefore, c ~ x, and (n) ~ y^. For the same reason, on approaching the extreme 

delocalized regime some power-law dependencies on c ~ a; are recovered, which are typical to canonical ensembles of 
"full" random periodic matrices: e.g., ~ x^^^, a^^^ ~ x"^/^ and Sn = const. 

Though the general approach described in Section II is still valid, some modifications of the picture described for 
the Elliptic ensembles are imposed by the specific features of the PBRM ensemble; in the first place, by the symmetry 
breaking occurring ai ip — 0, due to which this ensemble is not stationary. One-point spectral statistics (level spacing 
distributions, curvature distributions, ...) change as ip is increased from to some value ipc which is however small, 
in the sense that it decreases with 1 / y^; at large c the transition is therefore a very sharp one . The correlation 
analysis of Sec. II cannot be literally replicated for the PBRM model, as the distribution of BPs in the complex 
(/?— plane cannot be homogeneous in ip; the width 7 of a complex branch-point may not be statistically independent on 
its location ipo, because of the changes of symmetry occurring at special values of ip. We have, therefore, to investigate 
how the statistical distribution of BPs is affected by the symmetry breaking. This is in turn closely related to the 
distribution of avoided crossing sizes (gaps) and their positions, that will be analyzed below. 

B. Avoided Crossings. 

The above mentioned symmetry of the PBRM model entails that eigenvalues e{ip) are symmetric with respect to 
ip = 0,7T ((/3 = tt), and their velocities vanish there. It follows that relative extrema for the level spacings are observed 
at these symmetry-breaking values of ip. In other words, a large set of avoided crossings (relative minima of spacings) 
must appear at a fixed value of (/? = (or tt). Thus, with probability one, there will be BPs lying on the lines ip — 0,tt. 
Such BPs will contribute a delta- function term in the distribution of real parts ipo of BPs. Moreover, since the width 
of any of these avoided crossings is at once a nearest-level spacing, the distribution of these widths, at small values, 
will behave exactly like the level spacing distribution at </? = 0,7r, hence with the exponent (3=1 (and not (3 — 1, as 
with BPs having a random position). 

We therefore infer that the symmetry breaking lines contribute a singular term to the distribution ^(7, ip) of BPs 
in the complex plane, of the form (apart from a normalization constant): 
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(4.2) 



with /o(7) ^ 7 at small 7. We shall now examine the other type of avoided crossings, namely those which have 
fluctuating positions inside the (0, tt) interval. In order to numerically compute the positions and sizes of avoided 
crossings we have used the same procedure as in previous studies on stationary ensembles; we refer the interested 
reader to (24| for details. 

In fig.^ we illustrate the distribution of the positions of the "fluctuating" avoided crossings lying inside (0,7r), for 
matrices belonging to the diffusive regime. A strong dependence of the abundance of avoided crossings with respect 
to (f is apparent. In particular, fluctuating BPs appear to be repelled by the symmetry breaking values (where, as 
remarked above, other BPs already exist.) . The abundance of fluctuating avoided crossings depends on the position 
</? in a way which resembles that of the average squared level velocity, which will be discussed later. 

Although the distribution of avoided crossings is not homogeneous in ip, the (conditional) distribution of their sizes 
appears to be the same at all 7^ 0, tt and is presented as an insert in fig.^. Thus we argue that the BPs located inside 
(0, 7r) contribute a "smooth" part po{'y,ip) to the total density of BPs, which is not iy9-independent in the vicinity of 
0, TT, but is ^ 7 at small 7, at all ip. 



C. Spectral Correlations 



On the grounds of the above discussion, the distribution of the BPs in the complex plane consists of a " smooth" 
part, given by the BPs internal to (0, tt), plus a "singular" part, exactly located on the symmetry breaking lines. In 
both parts, the probability of small "f's scales with the same exponent 1. The correlation analysis of Section II now 
applies unaltered to the effect that the type of singularity exhibited by Cab(v) at (/j = is of GUE type, like that of 
(72(1^9); the numerical coefficients, nevertheless, may be different (and will indeed become different on moving towards 
the localized regime). 

On numerically computing Fourier transforms, we have in fact found distributions c„ whose general features cannot 
be distinguished from those observed in the pure GUE case, on the level of accuracy of our computations. In 
particular, a decay c„ ~ is apparent (figs. |l0|,|l2| ). The corresponding correlation functions arc almost identical to 
those obtained for the Elliptic GUE case (flg.(|ll|). 

The n""* decay is by no means consistent with the behaviour suggested by Braun and Montambaux pO| , which 
is instead of the type exhibited by Ci{ip) in ( 2.17 ) and would therefore result in a decay. Berry and Keating 
(BK) have proposed an explicit semiclassical approximation for Cab-, which is compared with the numerically 
computed CF for the PBRM model in flg.|l^. A free parameter, w* , which appears in the BK formula, has been fixed 
by means of a one-parameter fit to our data for N = 201, b = 64. The obtained value w* — 3.92 is quite small; 
therefore, we have used Eq.(28) of ref. ||l^, instead of its limiting form, Eq.(29), which is valid in the limit of large 
w* . Including the "non-diagonal" correction [ pT| yields no significant improvement to the fit. The minimum of the 
semiclassical CF is too shallow as compared to numerical data (see inset in fig.|To|). Moreover, since the Cab of BK is 
an analytic function at all real ip, itsJourier transform has an exponentially decaying tail, at variance with the real 
algebraic tail (main panel in fig.pi]|). 

It is worth pointing out, as a purely empirical observation, that a significant de facto agreement with numerical 
data is provided by a recently proposed semi-empirical representation of the CF p6[ | , which we show in the same 
inset. No fitting was used in this case, the only adjustement being the rescaling from the dimensionless parameter X 
(for which the semiempirical representation of CF was given) to (p. Different data sets (e.g. for N = 201, b = 101 , 
close to the ballistic regime, with w* ~ 6.35) have shown an even better agreement with this analytical representation 
(while no significant improvement has been observed with the BK formula). 

Nevertheless this semi-empirical representation , too, is in the form of an analytic function, and cannot therefore 
correctly represent the small- ip behaviour. Like the BK formula, in the Fourier domain it yields exponential instead 
of algebraic tails, at variance with numerical data. 

On moving towards the localized regime, the C((p) curve undergoes a continuous deformation (flg.p^, and the 
singularity at ip = Q remains of the same type all the way down to the deeply localized regime; in fact, Fourier 
transforms always exhibit a \n\~'^ tail. In all the investigated range we have obtained evidence that C{(p) depends 
on b,N only through the localization ratio x = b'^/N: C{(p,b,N) « C{(p,x) at b,N >> 1. In the diffusive regime, 



^The value of w* can also be estimated from the large (p behaviour, Eq.(39) of ref. a much larger value is found in this 
way. Though the large (p behaviour is now correct, the small and intermediate (p behaviour is in even stronger disagreement 
with numerical data. 
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x > 1, the scaling (^) is still approximately valid: C{ip,x) « c{x)^AB{Tr'P\/c{x)) with c{x) = C{0,x) x (weak- 
localization corrections are hardly assessed at the present level of accuracy). As an illustration, we plot in fig.jl^) 
the position (po{x) of t he ze ro of C{ip,x), as a function of x. In the diffusive part of the plot, foix) is proportional to 



^/x, consistently with (1.2) and with c ~ x. In the localized regime this scaling is not valid any more. The Fourier 



spectrum shrinks exponentially fast; all the harmonics except the 1st tend to disappear, so that C{ip) tends to ccos(p. 

D. Statistics of Fourier coefficients 

In the PBRM ensemble, the Fourier coefficients a„ are real. Our numerical data for their statistical distributions 
show in the delocalized regime a behaviour similar to the one observed for stationary ensembles. At small n < (n) 
we obtained a Gaussian distribution of a„; as a consequence, the log - variance af^ (i.e., the variance of In (a^)) is 
constant and equal to 7r^/2. . For n larger than (n) the distribution exhibits broad, exponential tails (fig.(|l4|)). 

On the other hand, in the localized regime, none of the coefficients is Gaussian-distributed. 

In the stationary case, the assumption that the phases ipoj of the BPs come at random immediately leads to the 
vanishing of correlations between different a„ - a result which we have obtained in a different way, and which is 
intimately linked to the stationarity of the ensemble. In the PBRM case, which is not stationary, different Fourier 



coefficients have nonzero correlations; the asymptotic formula ( |2.9| ) can be used together with formula (2.6) to study 
the correlation matrix (anOm)- 



The random phase assumption is not legitimate for the BPs which have nonrandom phases 0, tt. From (2.6),(2.£) 
it follows that these BP contribute a correlation 

-7(|n|-|-|m|) 1 

(a„a™) ~ (1 + (-l)"+™)(7^ ) ^ (1 + (-l)"+™)Af- (4.3) 

|nm|2 |nrn| 2 -I- |m|)'^ 

with M a numerical factor. In estimating the ensemble average, the GOE-like distribution of small 7 of BPs lying 
on the symmetry line s wa s used. We shall neglect the contributions of the remaining BPs, which have fluctuating 



phases. The formula (4.3) will be used in the next Section to analyze the behaviour of certain quantities near the 



symmetry-breaking point ip = 0. 

E. Kinetic Energies and Curvatures. 



From the expansion 



we obtain, using ( |4.3| ), 

k 

where the coefficients have the asymptotic behaviour: 



^ max{\k\,\p\Y \k\^ 

As we know, such an asymptotic behaviour of the Fourier expansion entails a singularity dX ip — 0, of type (e''^{(p)) ^ 
(/?^ln|(^| (and a similar behaviour close to (p = tt), as predicted in for a GOE-GUE transition and numerically 
observed for the Anderson model |10| ]30|] . 

A quite similar analysis shows, that the 2nd moment of the curvature has a logarithmic singularity at = 0, tt, as 
again it was predicted in [pl|,p^D 



*We note in passing that this singularity is an integrable one, and that t his fa rt alone is sufficient to exclude that the 2nd 
derivative of C{ip) may be singular at tp = 0: see the remarks following eqn.|2.20|. 
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Figure shows the dependence of the average squared velocity (normaUzed to its maximal value ) as a function of 
the scaling variable, X = ttlp^/c^ for N = 71 and for several band widths b. Our numerical data are consistent with 
the mentioned logarithmic behaviour. Moreover, in the delocalized regime, the data for different values of b follow 
the universal dependence recently predicted for the GOE-GUE transition |Q 

(e'((^)2) _ _^ yV2crfi(%/2X)[l + 2X2exp(2X2)Ei(-2X2)] 

, (4-4) 



e'la. 2Xexp{2X^ 

where Ei(.) denotes the exponential integral and erfi(x)=-i erf(ia;) is the imaginary error function ]29t . 

The statistics of the mean kinetic energy is similar to the one discussed for Elliptic ensembles; again, the distribution 
P{r) is better approximated by a xt distribution than by a log-normal one. In the localized regime, the width of the 
distribution P[ln(r)] increases with decreasing x. As shown in figj^, on decreasing x the logarithmic variance cr^i^ 
grows faster than x~^^'^. 

The relation ( 2.20| ) is still valid, as is the proportionality between the width of the curvature distribution P{K) 



at iy9 = 0, measured by gi ~ (|_ft"(0)|), and c [g7| The proportionality constant is close to 27r, as numerically 
observed for the Anderson model [0 and analytically found in j2^. Localization effects produce deviations from 
both proportionality relations at cc << 1 |^ 

Within the above described Fourier framework, it is also possible to analyze how the distribution of the curvature 
K(ip) depends on ip. This problem is also attracting interest; here we shall just make some qualitative remarks, 
deferring a more detailed analysis to future publications. The curvature distribution is affected by the distribution of 
BPs in the vicinity of the chosen point if on the real axis; in fact BPs coming very close to this point determine large 
curvatures, hence their statistical distribution determines the tail of the curvature distribution. However, the singular 
and the smooth parts of the BP distribution contribute in a different way in this tail. The singular part yields a 
relatively large probability for BPs coming close to ip = 0, tt, which in turn yields a large probability of large curvatures 
there; as a matter of fact, the curvature distribution at ip — O^n was found to be of the Zakrzewski-Delande type 
Nevertheless the same BPs, being bound to the lines (/? = 0, tt, cannot come arbitrarily close to a point 93 ^ 0, tt, 
and cannot therefore produce large curvatures there. In fact, if only the singular part of the BP distribution is taken 



into account, then a simple estimate based on (2.6),(|2.9|), (4.3) shows that at 7^ 0,7r all the moments of K{ip) are 
finite. Therefore, in any ensemble for which a smooth part is missing, the tail of the curvature distribution at 7^ 0, tt 
decreases faster than any power of curvature. This is exactly the case with the 2x2 matrix ensemble considered in 
pl| , where all BPs lie at ip — {j ~ in the notation of |^l|), and the curvature distribution was in fact shown to 
have Gaussian tails at (/? 7^ 0. 

Short of being typical of ensembles exhibiting a magnetic flux-induced GOE-GUE transition, such a behaviour rather 
appears as an artifact of ensembles lacking the smooth part of the BP distribution; the presence of such a part would 
in fact restore a chance for BPs coming close to any point on the real line. The impact of this fact on the curvature 
distribution is currently under investigation. 

V. CONCLUSIONS. 

In this paper we have illustrated the usefulness of analytic continuation to complex parameter values in the analysis 
of various statistics of parametric level dynamics, which are essentially determined by the statistical distribution of 
complex level crossings, ie, of the singularities of branch-point type exhibited by the eigenvalues as functions of the 
complex parameter. Though this approach is not, for the time being, as powerful as other methods (such as, e.g., 
supersymmetric ones), it offers, in our opinion, a particularly perspicuous picture. 

Our analysis was in many respects a heuristic one, although one rigorous result was also given: independently of 
symmetry, the level velocity correlation function can be analytic only if the matrix ensemble is so built, as to make 
too narrow avoided crossings impossible. Thus in most cases this correlation function has singularities, which have 
the same origin as the divergence of the 2nd moments of certain level derivatives, to which they are in fact connected 



by the relation ( 2.20 ). The latter is in a (rather superficial) sense a generalization of the Akkermans-Montambaux 
relation for conductance. 

Our estimates for the singularities of the correlation function were confirmed by numerical simulations, and also 
cover the physically important case, in which the parameter has the meaning of an Aharonov-Bohm flux. This case 
was modelled by the Periodic Band Random Matrix ensemble. Unlike most of the usually considered ensembles, 
which consist of "full" random matrices, this ensemble exhibits a proper diffusive regime in a suitable parameter 
range, marked by an ohmic dependence of the average curvature on the matrix size; therefore the corresponding 
eigenf unctions, although delocalized, have to differ from those of "full" random matrices in some as yet not fully 
understood, but nonetheless essential respect. The investigation of this ensemble has confirmed that various scaling 
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properties, originally established for conventional ensembles, remain substantially valid in the diffusive regime (already 
investigated in |l0| for a different model). 

The GOE-GUE breaking of symmetry which occurs on switching on the Aharonov-Bohm flux appears to have no 
impact on the behaviour of the correlation function, which looks the same as for a homogeneous GUE ensemble. 
Nevertheless, the transition is reflected by a singular behaviour of other quantities. Some features of this transition, 
such as e.g. the way it affects the curvature distribution, are not yet understood, as they may not be correctly 
reproduced by 2 x 2 matrix ensembles, which have a nongeneric distribution of branch points. 

For the case of stationary ensembles, we have shown numerically that, whereas the Fourier amplitudes a„ of the 
eigenvalues are approximately Gaussian distributed when n is small in comparison to the width of the Fourier spectrum, 
strong deviations from Gaussian distribution occur at large n. This behaviour can be qualitatively connected with 
the statistical distribution of level derivatives, which have Fourier expansions in which the role of the non-Gaussian 
amplitudes a„ is enhanced, the more, the higher the order of the derivative. While 1st order derivatives (level 
velocities) are still Gaussian distributed, this is no longer true already for 2nd order derivatives (level curvatures), 
which have in fact a Cauchy distribution; a strongly non-Gaussian distribution is therefore to be expected for all 
higher-order derivatives. 

Our theoretical apparatus heavily hinges on periodicity, and it is a natural question, whether any of our conclusions 
can be generalized to the general case of a non-periodic parameter dependence. Some kind of a stationary behaviour 
is needed, in order that a velocity correlation may be defined; an essential step is then to filter out any "secular" 
component of the eigenvalue motion, as can be done, e.g., by the unfolding process. In any case, the resulting branch- 
point pattern and the corresponding Fourier analysis may be, in our opinion, significantly different from the periodic 
case discussed in this paper. 
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VI. APPENDIX 

We compute analytically the velocity correlation for the level dynamics of 2x 2 matrices H{(p) = Hi cos(p + H2 sm(p, 
where 



■^[X2-iX3) a2 I \^{y2-m) 02 



both belonging to GUE or GOE ensembles; in the latter case X'i — yj, — 0. All matrix elements are normally 
distributed with the same variance. Setting ai + 02 = 's/Sa, 61 -I- 62 = %/2&, oi — a2 = ^/2xi and 61 — 62 = V^yi, we 
have the following expressions for the eigenvalues: 

£±(■^5) = —={acosip -\- bsaiip) ± —j=\xcofiip — ysimp\ (^-2) 
V 2 V 2 

We begin by computing the following eigenvalue correlator, which actually depends on the difference ip = Lp2 — 1^1: 

Piv) = {^Y^' j rfad6d"xd"ye--('^'+'''+l-l'+l^l')i[e+(^i)e+(^2)+e_(¥'i)e_(^2)] (A3) 

with n = 2, 3 respectively for GOE and GUE. The velocity correlator is then obtained as a second derivative: 
C(^) = -F"(v.)A-2. 

The eigenvalue correlator is the sum of two contributions Fi + i^2 7 where 

dadbe~'^'^°' ^ -iacosipi + 6 sin (/Si) (a cos ip2 + 6 sin 932) — -r- cos ip 
2 Alo 



F2 — {^^ j d^xd^ye "^'"^l^+l^l^^ilxcosyji — ?7siniy9i| • |xcosi^2 — y sin(p2| 
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With the linear change of variables, for different from and ±7r, 



{x cosipi — ysiuifi) , V = — {x cos (p2 — y sin ip2) 



sm ip sm if 

the second integral simplifies to 

F2 = -^|sin(pr+2 / (PuM I rf"^;|^;|e-l"l'-l''l'+23-^7cosv. M 4) 

27r"u; J J 

Note that F2 is an even function of if. 

Let us first consider the case of GOE. In polar coordinates, letting 6 be the angle between the two vectors u and v, 
we have 

F2{ip) = 2^ sin<^42^ J- u^duj^ v'^dve-^^-^" /o" d^e^"" ^ cos e 
= ^sin''(y9 J^u'^du j'^ v'^dve~^ ~'" Io{2uvcos(p) 

where Io{x) is a Bessel function. The double integral is tabulated [Prudnikov, vol. 2]and we obtain: 

F2('p) = — \2E(cosLp) - sin^ (fK(cos(p)] 

Alo 

where E and K are the complete elliptic functions. Adding the Fi term, and computing the second derivative we get 
the velocity correlator for GOE. Explicitly, it is a function oft = cos(p: 

^ = ^[t^+E{t)-{l-t')m] (A5) 

The small <^ expansion has the expected logarithmic term: 

^('^) 1 _ ^2 ^^og 2 - i - ^ log 1^1) + 0{^') iA.6) 



C(0) ^ V2 8 4 

At (f = tt/2 {t = 0), the normalized correlator has the value — tt/S, and at = ti" — e it vanishes as — e^[(l/8) + 
(3/4)log(4/e)]. 

Let us now consider the GUE case. In spherical coordinates, letting 9 be the angle between the vectors u and v, 
we have: 

F2{'fi) = sin(/j|587r2 /~ u^duj^ v^dve-""^-"^ sin ^e^"'"^"^ cos e 

= r dujr dvV^e—^smH2V^cos^) ^^'^^ 

The double integral can be related to the following, available in the tables [Prudnikov, vol 1 pag 573] {pq > 1/4): 

f dx / dycosh(y^)e 
Jo 

For < (f <n the result is 



1" 7 73-7:7 arcsm , 

JO 4m -1 (4pg- 1)3/2 V2VM 



= 2^; 



/TT \ / sin^wN „ . 
7: — <^ 3 cos u) H + 3 sin M 

V2 V V cosy y ^ 



After summation of the Fi term, and derivating twice we obtain the velocity correlation for 2 x 2 GUE matrices: 



C(0) V2 Try TT Vcos</9 cos^t^y \2 2 

In particular, we compute the following Mac Laurin expansion: 



C(0) 37r' 

At = 7r/2 the function takes the value — 4/(37r) and vanishes as — (3/2)e^ iov ^p = -k — e. 
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Figure Captions 

FIG. 1. Multiple complex path used in the calculation of the Fourier amplitudes a„ with n > 0. Three branch points are 
shown. 

FIG. 2. Illustrating the scaling law (^^) in the Fourier domain for the Elliptic GUE ensemble. The universal function P2{t) 
[ Eq.(2.14)] versus t = n/{\n\) for matrices of rank TV = 71 (□), iV = 300(A), iV = 400(0)- Lines show the initial linear rise 
and the tail. 

FIG. 3. Large-n decay of the Fourier coefficients of the Correlation Function, for the three Elliptic ensembles: GOE(O), 
GUE(o), GSE(A). The Fourier coefficients have been divided by a/(C(0). 

FIG. 4. Scaling functions ^^{X) versus the rescaled parameter X for Elliptic ensembles and different matrix rank A'^: GSE 
with iV = 80 (full triangles); GSE with iV = 90 (•) ; GUE with iV = 71 (A); GUE with = 300 (Q); GUE with N = 400 
(□)• 



FIG. 5. The Scaling function ^2{X) for the GUE Elliptic ensemble, matrix rank 71. The dashed line is the 2nd order 
expansion 1 — 2X^; the solid line includes a 3d order term +1.033X"^. The insert illustrates the X~^ behaviour at large X. 

FIG. 6. Variance af^ of ln(|a„p) versus of t = n/{n) for the GUE Elliptic ensemble and iV = 40(»), 60(a) and 80(A). The 
horizontal dashed line corresponds to (complex) Gaussian distribution of an 

FIG. 7. Distribution P(lnr) of the logarithm of the mean square current, for the GUE Elliptic ensemble, obtained from 200 
realizations of rank A'^ = 40. The dashed line represents the log-normal distribution with {ln(r)) = 1.17 and crf,-j — 0.075; the 
solid line corresponds to the xl distribution, with mean c = 3.34 and number of degrees of freedom u ^ 

FIG. 8. Illustrating the power -law dependence of various ensemble averages on the scaling parameter x = /N , for Periodic 
Banded Random Matrices of rank 71,101,201. Average squared current c = (r) (o), mean current variance (□), log-variance 
""rin (*) 1 width of the Fourier spectrum (|n|) (A). Lines represent the limit dependences that are approached in the delocalized 
regime. The vertical line is somewhat arbitrarily drawn to represent the transition from localization to delocalization. 

FIG. 9. Abundance of avoided crossings as a function of for the PERM ensemble ( 150 matrices with A'^ = 201, 6 = 64). The 
histogram has been obtained from 4452 avoided crossings occurring inside the open (0, tt) interval , in which only eigenvalues 
belonging in the central half of the spectrum were involved. The dashed line is proportional to the average squared level velocity. 
The inset compares the distribution of avoided crossing sizes (normalized to a unit mean gap) with the Wigner distribution. 

FIG. 10. CF for the PERM ensemble in the diffusive regime {N = 201,6 = 64) (inset, full line) and its Fourier transform 
(main panel, full circles ). The dotted line represents the decay. Both in the main figure and in the inset, dashed lines 
represent the Berry -Keating semiclassical approximation. The dotted line in the inset represents a semi-empirical formula 
proposed in ref. [jl6|. 



FIG. 11. Comparison of the scaled correlation functions $2 for the Elliptic GUE ensemble with A'' = 71 (circles) and <^aB 
for the PBRM emsemble with N = 101, & = 51 (shown by a continuous line, obtained by interpolating numerical data). 

FIG. 12. the Correlation Function for the PBRM ensemble (a), and its Fourier transform (b), for various values of the 
localization parameter. Numerical data in (b) are for matrices of rank A'^ = 201 and b — 101, 64, 23, 9 from rightmost to 
lefmost curve. All the reported data sets have the same slope, even in the localized regime (for fe = 9, \n{b'^/N) -1). 
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FIG. 13. The position (po of the zero of the CF for the PERM ensemble, versus the localization parameter x = /N for 
various matrix sizes. The dotted curve corresponds to ifio = 6.55x^^^ /n. 

FIG. 14. Semilogarithmic plot of the distribution of the coefficient n^ttn for the PERM ensemble, with N = 61,6 = 31, and 
forn = 1 (□), and n = 10 (•). 

FIG. 15. Average squared velocity — {e'{(p)^) (normalized to its maximum value) versus the rescaled phase X , for the 
PERM ensemble with N = 201 and b = 9(A), 23(o), 64(D) and 101(»). In the delocalized regime numerical data agree with 
the theoretical formula of Taniguchi el al, (solid line). The insert shows the same data without rescaling. 
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